Manipulation of the dynamics of many-body systems via quantum control methods 
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We investigate how dynamical decoupling methods may be used to manipulate the time evolution 
of quantum many-body systems. These methods consist of sequences of external control operations 
designed to induce a desired dynamics. The systems considered for the analysis are one-dimensional 
spin-1/2 models, which, according to the parameters of the Hamiltonian, may be in the integrable 
or non-integrable limits, and in the gapped or gapless phases. We show that an appropriate control 
sequence may lead a chaotic chain to evolve as an integrable chain and a system in the gapless phase 
to behave as a system in the gapped phase. A key ingredient for the control schemes developed here 
is the possibility to use, in the same sequence, different time intervals between control operations. 



I. INTRODUCTION 



' The study of one-dimensional quantum many-body systems out of equilibrium offers several fundamental challenges. 
A main question is whether an isolated quantum system can or not thermalize after a quench. While a relationship 
between quantum chaos and thermalization is certain to exist [11-0], an ongoing controversy is the occurrence of 
thermalization for a chaotic system in the gapped phase In terms of transport properties, the discussions are 

associated with the expected dynamical behavior of systems in the integrable and nonintegrable regimes [Tol - [20| , and 
in the gapped and gapless phases [2ll,[23|. Extraordinary technological advances have motivated the use of experiments 
^ to address these issues. In this work we show how these experiments may benefit from dynamical decoupling methods 

to accurately reach the different regimes and phases of interest. 
^ The fundamental questions raised above have found a test-bed in experiments with ultracold quantum gases in 

O^i optical lattices 23] and with magnetic compounds [24- 26]. Both systems require a high level of controllability. In 
optical lattices, control of atomic transport is achieved with time-dependent variations of the depth and phase of the 
y—{ . lattice potential; whereas in magnetic systems the preparation of samples where impurities and unwanted couplings 
^ ' are negligible is key for obtaining a desired transport behavior. In addition to monitoring experimental conditions, 
CO , we advocate the use of sequences of external electromagnetic field pulses to manipulate the dynamics of the systems. 
■ This method of quantum control has long been employed in nuclear magnetic resonance (NMR) spectroscopy |27l [28| . 
\ There, the Hamiltonian of nuclear spin systems is modified with sequences of unitary transformations (radio-frequency 
. ' pulses), which constantly rotate the spins and remove or rescale selected terms of the original Hamiltonian. These 
CO ' techniques have acquired a new dimension in the context of quantum information processing |2& . i30i] , where the studies 
of transport of information [31] and the manipulation of electron spins with ultrafast laser pulses [32, [33 are among 
the latest developments. Recently, these methods have been proposed also for the control of atomic transport in 
spinor optical lattices [3^ . 

In the present work we use the theory developed in NMR to construct sequences of control operations that appro- 
priately modify the Hamiltonian of a one-dimensional quantum many-body system and forces it to follow a desired 
; I ■ dynamics. The systems analyzed are modeled by spin-1/2 chains, which may also be mapped onto systems of spinless 
d ' fermions or hardcore bosons. We focus on limits that are of interest for studies of transport behavior and thermal- 
ization. Specifically, our goals are (i) to dynamically remove the effects of terms leading to the onset of chaos, such 
as disorder and frustration, and compel the chain to evolve as an integrable system, as introduced in Refs. [l8lll9|: 
(ii) to dynamically increase the effects of terms that can drive the system across a transition from the gapless to the 
gapped phase (and vice- versa), aiming, for example, at forcing a conductor to behave as an insulator or at bringing 
a chaotic system into the gapped phase. For the gapless-gapped transition, the strengths of the coupling terms in 
the Hamiltonian are modified by applying to the target system sequences with different time intervals between the 
control operations. The idea of varying the time separation between pulses has been exploited in sequences created to 
suppress the couplings between a qubit and its environment [H,!!!]. Contrary to that, the systems considered here 
are isolated from the environment and our objective is not necessarily to completely eliminate couplings, but instead 
to control their strengths. 
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This paper is organized as follows. Sec. II describes the model under investigation. Sec. Ill reviews a technique to 
design control sequences, which is based on the average Hamiltonian theory. Sec. IV presents in detail sequences that 
dynamically lead the system studied across a transition between the nonintegrable and integrable limits and between 
the gapless and gapped phases. Numerical results are shown for the time evolution of local magnetization based on 
a particular initial state. A discussion about control metrics and the dependence of the outcomes on initial states is 
left for the Appendix. Conclusions are given in Sec. V. 

II. DESCRIPTION OF THE MODEL 

We study a one-dimensional spin- 1/2 system with open boundary conditions described by the Hamiltonian 

Ho ^ Hz + H^N + chHt^-nn, (1) 

where 

L 

Hz — ^nS^, 
n=l 
L-1 



n=l 
L-2 



Above, ft is set equal to 1, L is the number of sites, and Sf^'^'^ = a^'^'^/2 are the spin operators at site n, cr^'^'^ being 
the Pauli matrices. 

(i) The parameter e„ corresponds to the Zeeman splitting of spin n, as determined by a static magnetic field in the 
z direction. The system is clean when all e„'s are equal (e„ = e) and it is disordered when at least one spin n has an 
energy splitting different from the others (e„ ^ e). We therefore consider only on-site disorder. 

(ii) -ffNN, also known as the XXZ Hamiltonian, corresponds to the nearest-neighbor (NN) exchange; S^S^^i is 
the Ising interaction and S^S^j^^ + S-^S^j^^ is the flip-flop term. The coupling strength J and the anisotropy A are 
assumed positive, thus favoring antiferromagnetic order. By varying A in Hq = Hjqj^ we may go from the gapless 
phase (A < 1) to the gapped phase (A > 1), which, at T = 0, is followed by the metal-insulator Mott transition 
[331 . In the particular case of A = 0, we deal with the XY model; whereas A = 1 leads to the isotropic Heisenberg 
Hamiltonian (38| . 

(iii) The parameter a refers to the ratio between the next-nearest-neighbor (NNN) exchange, as determined by 
-f^NNN, and the NN couplings. The inclusion of NNN antiferromagnetic exchange frustrates the chain, since NN 
exchange favors ferromagnetic alignment between the second neighbors. Different phenomena are expected in the 
presence of frustration |39f]. In the case of a closed isotropic system, for instance, a critical value ac ~ 0.241 exists 
which separates the spin fluid phase (a < ac) from the dimer phase (a > Uc) [40ll4l|. 

TABLE I: Parameters of the spin- 1/2 chain described by Eq. ([l| 



£71 


Energy sphtting of spin n 


J 


Coupling strength 


A 


Ratio between Ising interaction and flip-fiop term 


a 


Ratio between NNN and NN couplings 



A clean one-dimensional spin-1/2 model with NN exchange only is integrable and solved with the Bethe Ansatz 



method |42|. The addition of defects [43|, even if just one in the middle of the chain [4J|, or the inclusion of NNN 
couplings |45l. l46j may lead to the onset of quantum chaos. 

III. DESIGN OF DYNAMICAL DECOUPLING SEQUENCES 



Dynamical decoupling methods consist in the application of sequences of unitary transformations to generate a 
desired effective propagator at time t. The first experimental implementations of sequences of electromagnetic pulses 
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aiming at controlling the dynamics of molecular systems appeared in the context of NMR spectroscopy and correspond 
to the so-called spin-echo and Carr-Purcell sequences [U I3 . The basic idea of the method is to add a time-dependent 
control Hamiltonian, Hc{t), to the original (time- independent, in our case) Hamiltonian of the system, Hq, so that 
the system evolution becomes dictated by the following propagator 

U{t) ^TexpS^-i [Ho + H^{u)]duY 

where T stands for time ordering. This operator may also be written as [13) HI], 

U{t) = Uc(t)Tcxp |-i ^ [U-\u)HoUc(u)]duj , 



where 



Uc{t) = Texp 



-i / Hc(u)d' 



Jo 



depends only on the introduced perturbation Hdt) and is known as the control propagator. We deal with cyclic 
control sequences with cycle time Tc. This means that the perturbation and the control propagator are periodic, that 
is, Hc{t + nTc) — Hc{t) and Uc{t + nTc) — Uc{t), n G N. Since Uc{0) — 1, it follows that the control propagator at 
T„ = nTc is Uc{nTc) = 1. As a resuh, the general propagator at Tc is simply U{Tc) = re"* /o^'['^e"'(")^o'^-(")l'*" and 
for n cycles 

U{nTc) = UiTcT. 

Therefore, the evolution of the system at any time r„ — nTc requires only the knowledge of the propagator after one 
cycle [23,|2i|. 

We consider the ideal scenario of arbitrarily strong and instantaneous control operations, known in NMR as hard 
pulses and popularized in the field of quantum information as bang-bang control (29j . Each pulse Pk+i is applied 
after a time interval of free evolution t^+i = ife+i — ifc, where A: g N and tg = 0. For a cyclic sequence with to pulses, 
the evolution operator at Tc = nir is then given by, 

)...P2U{t2,ti)PiU{ti,0) 

= {PrnPrn-1 ■ ■ ■ A) (-P,n-1 ■ ■ ■ Pl)^ U {tm , tr,r-l){Prn-l ...Pi) .. ■ (P2P1 ) ^ t/ftg , ^2 ) (P2P1 ) £7(^2 , ti )PiJ /(^ i , 0) 

Uc{Tc) = 1 g-i(P„-l...Pl)tffo(Pm-l...Pl)T„ g-i(P2Pl)*ff0(-P2Pl)T3 g " i P/ ^0 Pi T2 g-iHoTl 

= exp [-ii?„i_iT„i] . . . exp [~iH2T3] exp [-iHiT2] exp [-iHaTi] 

= exp [^iHTc] . (2) 

Above, the notation Hm~i ~ (Pm-i • • ■ Pi)^-ffo(Pm-i ■ • • Pi) has been used. The Hamiltonian H in the last line is 
obtained via Baker-Campbell-HausdorfF expansion and is referred to as the average Hamiltonian p7l. [28j. 
The lowest order term of iJ in r corresponds to 



m— 1 
fc=0 

In designing a dynamical decoupling sequence, the primary goal is to guarantee the proximity of H^'^^ to the desired 
Hamiltonian. This is the focus of the current work. We note, however, that various strategies exist to eliminate some 
of the higher order terms [27l [28| or to reduce the accumulation of errors arising from imperfect averaging pol - ls^ . 



IV. MODIFYING THE SYSTEM EVOLUTION 



In this section we explore how dynamical decoupling schemes may induce transitions between different regimes and 
phases of a quantum many-body system. We aim at dynamically reshaping the original Hamiltonian ([1]) to make it 
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as close as possible to a desired Hamiltonian (the subscript "w" is used whenever we refer to parameters and 
quantities associated with the wanted evolution). Our numerical results are performend in the full Hilbert space of 
dimension D — 2^ . We cannot take advantage of the symmetries of i/o llj, such as conservation of total spin in the z 
direction, to reduce this dimension, because the intervals of free evolution of the pulsed system involve Hamiltonians 
\Hk in Eq. ([3])] that are different from the original Hq. 



A. Non-integrable vs integrable systems 



To start, we discuss sequences of control pulses that average out the terms in Hamiltonian ([IJ that may lead to the 
onset of chaos, namely on-site disorder and NNN exchange. 



1. Annihilating the effects of on-site disorder: + Hnn — > Hnn 

The first and simplest sequence that we present was proposed in [isl . and it aims at eliminating the effects of 
on-site disorder. Consider a disordered system described by Hz + H^n with one or more sites having e„ 7^ e. In order 
to recover the transport behavior of a clean chain given by Hnn, we apply, after every r, control operations P, which 
rotate all spins by 180° (7r-pulses) around a direction perpendicular to z. Around x, for instance, the pulses are 



Px = exp ^-i-TT ^ = exp(-i7rS'^). 



(4) 



The propagator at cycle time Tc = 2t becomes 

U{T,) = exp [-^ {H, + Hnn) t] P, exp [-t [H, + Hnn) t] 
= (-1) exp [-i {-Hz + Hnn) t] exp [-« {H^ + Hnn) t] , 

which leads to the following dominant term of the average Hamiltonian 

= Hnn, 

as desired. 

2. Annihilating the effects of NNN exchange: Hnn + Hnnn Hnn 

A sequence to eliminate the effects of NNN exchange (and also on-site disorder) was introduced in [3] and is further 
investigated here. It has eight 7r-pulses equally separated (Tc = 8r) and applied to selected spins, which therefore 
assumes site addressability. The pulses are 

L(L-1)/4J L(i-2)/4J 

fe=0 fe=0 
L(L-3)/4J L(i-4)/4j 

fe=0 k=Q 
L(L-2)/4J L(i-3)/4j 

P5 = Pr^ n e~"^2"+« Yl e~'''^3+4f= ^ 

fe=0 k=Q 
L(L-1)/4J L(i-4)/4j 

Pe = Ps= n e~"^i+« Yl e-'^^"+^K (5) 
fe=o k=a 

They change the signs of the coupling terms in each interval of free evolution according to Table |TT1 
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TABLE II: Sign changes of the coupling terms in Hnn and Hnnn according to the pulses ((5|. 



NN couplings 


(0,r) 


(r, 2r) 


(2r, 3r) 


(3r, 4r) 


(4r, 5r) 


(5r, 6r) 


(6r, 7r) 


(7r, 8r) 


rtx QX 
•^2/0 + 1 '-'2fe + 2 


+ 


+ 


+ 


+ 


+ 


+ 






QX QX 

'-'2fc+2^2fc+3 


+ 


+ 






+ 


+ 


+ 


+ 


qV qV 
2fc+l 2fe+2 


+ 


+ 


+ 


+ 


+ 


— 


— 


+ 


qy qy 


+ 






+ 


+ 


+ 


+ 


+ 


QZ QZ 

'~'2k+l02k+2 


+ 


+ 


+ 


+ 


+ 




+ 




QZ QZ 

'-'2k+2^2k+3 


+ 




+ 




+ 


+ 


+ 


+ 


NNN couplings 


(0,r) 


(r, 2r) 


(2r,3r) 


(3T,4r) 


(4r, 5r) 


(5r, 6r) 


(6r, 7r) 


(7r,8r) 


qx qx 
'5fc + l'-'fc+3 


+ 


+ 






+ 


+ 






qy qy 
"fc+l"fc+3 


+ 






+ 


+ 






+ 


Sk+lSk+3 


+ 




+ 




+ 




+ 





The sequence removes NNN couplings, but at the price of also reducing the effects of NN couplings. We achieve 
the following dominant term for the average Hamiltonian, 

Hio) ^ Hnn 
2 

In Fig. [1] we show the evolution of a chaotic chain described by Hn n + Hn nn under free evolution and subjected 
to the pulses ([S]), and compare it with the dynamics of an ideal integrable system given by Hnn/'2- The quantity 
considered is the local magnetization, which is defined as the magnetization of the first half of the chain, 

L/2 

M{t)^{^{t)\Y,S:,\^{t)). (6) 

n=l 

We choose as initial state |^'(0)) a highly excited state far from equilibrium: it consists of spins pointing up in the 
first half of the chain and pointing down in the other half, so that M(0) — L/4. A discussion about control metrics 
and the dependence of the results on different initial states is given in the Appendix. 




2 4 6 8 10 12 14 

Jt 



FIG. 1: Time evolution of local magnetization for a clean spin-1/2 chain with A = J/2 and L = 10. Blue dashed line: 
chaotic system -ffivjv + Hnnn in the absence of pulses; black solid line: integrable system Hnn/'2', symbols: chaotic system 
Hnn + Hnnn in the presence of pulses ((5|; triangles: r = 0.025J~^, squares: r — 0.05J~^, pluses: r — 0.1J~^. Main panel: 
data acquired every T„ = nTc = 8r, inset: data acquired every t„ = nr. 

In order to average out unwanted terms of the Hamiltonian, we need to consider cycle times smaller than the time 
scale of the free evolution of the system, that is, JTc < 1 ^5j|. For the T^s considered in Fig. [I] the agreement between 
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the evolutions of the chaotic system under control pulses and the integrable chain is excellent. As the values of r 
increase, higher order terms in the average Hamiltonian become significant, specially at long times, and the pulsed 
system must eventually separate from the ideal scenario. Even though strategies exist to reduce the effects of higher 
order terms in H j49l - l52| . the access to short time intervals between pulses is essential to the success of dynamical 
decoupling methods. 

Note that control sequences are designed to match the aspired evolution at the completion of cycles, but not 
necessarily in between cycles. This is well illustrated in the inset. The symbols indicate data obtained at every 
T„ = nTc, while the solid lines correspond to the values of M{t) at every tn = nr. The largest mismatches between 
the evolution determined by iJ^vAr + ^^iVATiv + Hc{t) and by occur in between cycles. 

A second alternative to recover the dynamical behavior of an integrable system consists in canceling out the NN 
exchange leaving the NNN couplings unaffected. This is achieved with a sequence of four pulses equally separated 
(Tc = 4r), which alternates rotations around perpendicular directions, say x and y. The pulses can affect only spins 
in odd sites, or only in even sites, or they alternate between odd and even sites as below, 

L(L-1)/2J L(i-2)/2j 

Pi=P3= Y[ e-^'^i+^S P2 = Pi= n e-*"^^^^'^. (7) 

fe=0 fe=0 

The pulses ([7]) change the signs of the terms of H^j^ + Hj^^j^ as shown in Table Hill 



TABLE III: Sign changes of the coupling terms in Hnn and Hnnn according to the pulses ([7]). 



NN couplings 


(0,r) 


(r,2r) 


(2r, 3r) 


(3r,4r) 


QX QX 


+ 






+ 


QV QV 

'-'fc+l'^fc+2 


+ 


+ 






QZ QZ 


+ 




+ 




NNN couplings 


(0,r) 


(r,2r) 


(2r, 3r) 


(3r,4r) 


QX QX 

'-'k + l'-'k + 'i 


+ 


+ 


+ 


+ 


QV QV 
•^k + l'^k + 'i 


+ 


+ 


+ 


+ 


QZ QZ 


+ 


+ 


+ 


+ 



The above control sequence leads to two uncoupled chains, one consisting of exchange couplings between spins in 
the original odd sites and the other one consisting of couplings between spins in the original even sites. The first order 
term in the average Hamiltonian is simply, 

= Hnnn- 



B. Gapped vs gapless systems 

Here, we analyze how the time evolution typical of a system with an energy gap may be achieved from a gapless 
system subjected to a quantum control scheme. Instead of completely averaging out unwanted terms of the Hamilto- 
nian, as done in the previous section, the goal now is to change the ratio between different coupling strengths. This is 
accomplished by varying the time interval between pulses, so that the weight of some specific terms may be decreased 
relatively to others. 



1. Increasing the role of the Ising interaction: A < 1 — !> A > 1 



Consider a given system in the gapless phase described by Hnjsi, where A < 1. To induce the transport behavior 
of a chain in the gapped phase, we develop a sequence that leaves the Ising interaction unaffected, but reduces the 
strength of the flip-flop term. The sequence consists of two 7r-rotations around the z axis, which affect only spins in 
even sites (or only in odd sites), that is. 



Pi = P2 = exp -ZTT V Sl,+, . (8) 
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The first pulse is applied after ri and the second after r2, which results in the following reshaped Hamiltonian at 

Tc = Tl+ T2, 

The gapped phase is obtained if TcA/{ti ~ T2) > 1. 

On the left panels of Fig. [2] we compare the desired evolution of a gapped system described by Hmn/'^, where 
= 2, with the dynamics of a given gapless system described by H^n^ where A = 1/2, and subjected to control 
operations. The pulses, given by Eq. have time intervals that lead to Tc/S./{ti — T2) = 2. The agreement between 
the two cases at the completion of each cycle is very good for the chosen time intervals, but of course it worsens as 
Tc increases. The price to dynamically achieve a transition from a gapless to a gapped Hamiltonian, however, is the 
slowing down of the system evolution in the gapped phase by a factor of 4. For a discussion about ballistic vs diffusive 
transport in XXZ models with small and large values of A, see Refs. [13, [13 ^.nd references therein. 



A /r 2.51 
^2.4 

2.3 
2.2 
2.1 
2 

I 

• \ 

1 - \ 
- 
-1 - 

FIG. 2: Evolution of local magnetization for a clean spin-1/2 chain with L = 10, described by Hnn with A — 1/2 (left panels) 
and by Hnn + aHNNN with A = 1 and a = 0.1 (right panels). Absence of pulses: blue dashed line; presence of pulses: symbols. 
Left panels: pulses ((8| leading to TcA/{ti — T2) ~ 2. Triangles (appearing only in the inset): n = 0.1J~^ and T2 = 0.06J~^; 
squares: ri = 0.2J~^ and T2 = 0.12J^^; pluses: ri = 0.4J~^ and T2 = 0.24J~^. Black solid line: system described by -ff]vjv/4 
with Aw — 2. Right panels: pulses Q leading to aTc/{T\ — To) = 0.64, where Ta — T2,3a ~ t-i(0.64 — Qf)/(0.64 + 3a). Triangles 
(appearing only in the inset): ri — 0.05J~^; squares: ri = 0.1J~^; pluses: ri = 0.2J~^. Black solid line: system described by 
{Hnn + Qiu-ffjv]vjv)/6.4 with A = 1 and — 0.64. Bottom panels: data acquired every cycle; top panels: data acquired every 
pulse. 

It is also possible to force a transition from a gapped to a gapless phase. For example, pulses (O separated by 
different intervals ti, T2 = T4, and T3, reshape the Hamiltonian of a Mott insulator described by HjsfN with A > 1 as 
follows, 

L-l 

H^'^ - E T [(^1 - {S:s:^, + SlSl^,) + A(ri + r3 - 2r2)5^5^+ J . 

The gapless phase is then obtained if A(ti + T3 — 2x2) < (ti — t^). In particular, when ti + = 2r2 we recover the 
XY model. 



2. From the spin fluid to the dimer phase: a < ^ a > Ua 

Let us now focus on a system described by H^ini + oiHninin with a < Uc and A = 1. Our goal is to apply to it 
a sequence of control pulses that induces a dynamical behavior similar to that of a system in the dimerized phase 



(n - T2) 
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(a > ac). A possible sequence consists of the four pulses presented in Eq. d?]), although now separated by different 
time intervals ti and T2 — = T4 = Ta- It leads to the following dominant term of the average Hamiltonian at 

= n + T2 + Ta + T4, 

^ c 

The desired dynamics is reached when aTd {ti — To) > occ- 

On the right panels of Fig. [5J we compare the dynamics of a system in the gapless phase given by Hj^pf + aff^vTVAfj 
where a = 0.1 and A = 1, with the wanted dynamics of a system in the gapped phase described by (-ffjvAr + 
aiJArjvAr)/6.4, whcrc = 0.64 and A = 1. The pulses ([7]) applied to the original system are separated by intervals 
that satisfy the constraint aTc/{Ti — Ta) — 0.64. For JTc < 1, the dynamics of the perturbed gapless system at the 
completion of the cycles agrees well with the evolution of the gapped system. 

Notice that our pulsed system will now be in the limit of strong frustration and will therefore become chaotic [46| . 
We expect an initial state with energy far from the edges of the spectrum to show a dynamical behavior typical of a 
chaotic system [53], although slowed down by a factor a^ja. 

A variety of new pulse sequences may be designed according to our needs and experimental capabilities. For 
instance, the sequence described above, which has no effect on NNN couplings, could be combined with a sequence of 
TT-pulses in the z direction given by 

L(L-1)/4J L(i-2)/4J 
fc=0 fc=0 

These pulses would affect only the flip-flop term of Hnnn- The combined sequence would reduce the strength of 
-ffjVAT and keep only the Ising interaction from H^nn. Hamiltonians similar to this resulting one have been analyzed 
in the context of transport [55] and thermalization in gapped systems Q. 



V. CONCLUSION 



We studied how quantum control methods may be used to manipulate the dynamics of quantum many-body 
systems. This is accomplished with sequences of unitary transformations designed to achieve a desired evolution 
operator. The unitary transformations correspond to instantaneous electromagnetic field pulses, which, in the case of 
spin systems, rotate the spins by a certain angle. The strategy employed to construct the control sequences was the 
average Hamiltonian theory. 

Four target spin-1/2 chains were considered. Two were chaotic: (i) disordered with NN exchange, (ii) clean with 
NN exchange and NNN frustrating couplings; and two were in the gapless phase: (iii) clean with NN couplings and 
small Ising interaction, (iv) clean with NN couplings and weak frustration. We discussed sequences of control pulses 
aiming at averaging out the effects of disorder, in case (i), and the effects of NNN couplings, in case (ii), and then 
induce in these chaotic systems a dynamical behavior typical of an integrable clean chain with only NN exchange. We 
also developed new sequences capable of changing the ratio between the coupling strengths of the system. This was 
done by exploiting, in the same sequence, different time intervals between the pulses, so that the effects of specific 
terms in the Hamiltonian may be reduced, but not necessarily eliminated. By dynamically decreasing the weight of 
the flip-fiop term in (iii), we obtained a time evolution characteristic of gapped systems with large Ising interaction 
(also mentioned was a scheme to achieve the dynamics of a metal out of an insulator). By decreasing the strength of 
the NN couplings with respect to the NNN exchange in (iv), we obtained the dynamical behavior of a chain in the 
limit of strong frustration (also mentioned was a way to combine to this scheme a sequence that removes the NNN 
fiip-fiop term). 

Site addressability and the possibility to vary the time intervals between pulses add versatility to dynamical de- 
coupling schemes and open an array of new opportunities for controlling quantum systems. Our focus here was to 
use these ideas as an additional tool to reach regimes and phases that are relevant in studies of the nonequilibrium 
dynamics of many-body systems. 



Acknowledgments 



J.D. thanks financial aid from the S. Daniel Abraham Honors Program at Stern College for Women. This work was 
supported by a grant from the Research Corporation. 



9 



Appendix A: Control Metric 



In the main text, we showed numerical results for the local magnetization. This choice was motivated by works 
on transport properties and thermalization, where quantities of experimental interest are usually considered. The 
decision to take the particular initial state where all up-spins are in the first half of the chain was a way to connect our 
results with previous transport studies p7l - [T9| . and also to provide a good illustration, since states where M(0) ^ 
would not give much information. However, in order to quantify how successful the designed control sequences are in 
general, we need an appropriate control metric which is independent of the initial state. 

How close a pulse sequence brings the evolution of a system to a desired dynamics given by Uw{t) may be assessed 
by the so-called propagator fidelity, which is defined as [56| , 



F^{t)^\Tr[Ulmim/D. (Al) 

Above, D = 2^ is the dimension of the system and U (t) dictates the evolution of the pulsed system. This quantity 
deals directly with propagators and therefore avoid any mention to particular initial states. Our goal is to achieve 
Fu{t) ^ 1. 




FIG. 3: Propagator fidelity vs time for clean spin-1/2 chains, L = 10. Left panel: propagator for — Hnn/^ vs propagator 
for a pulsed system described by -ffivjv + Hnnn', A = 1/2 for both. From top to bottom: pulses ((S| separated by r = 
0.025 , 0.05 , 0.075 , 0.1J~^ . Data acquired after every pulse, cycles indicated with symbols. Middle panel: propagator 
for i/m = Hnn/4: with A = 2 vs propagator for pulsed Hnn with A = 1/2. Data acquired after every pulse and indicated 
with symbols. From top to bottom: pulses ((HJ separated by ri — 0.1J~^ and T2 = 0.06J~^; ri — 0.2J~^ and T2 = 0.12J~^; 
n = 0.3J~^ and T2 = 0.18J~^; n = 0.4J~^ and T2 = 0.24J~^. Right panel: propagator for = [Hnn + a™-ff]viViv)/6.4 with 
aw ~ 0.64 vs propagator for pulsed Hnn + oHnnn with a — 0.1; A = 1 for both. From top to bottom: pulses ((Tjl separated 
by Ti = 0.025J"\0.05J"\0.075J~\0.1J"^ and r„ = T2,3,4 = ri(0.64 - q)/(0.64 + 3a). Data acquired after every cycle. All 
panels: fast decaying dashed line - absence of pulses. 

In Fig. 131 we show again the cases studied in Figs. [1] and [21 but now from the perspective of the evolution oper- 
ator. The performance of the pulse sequences deteriorates for longer T^s and longer times. For the shortest cycle 
times considered, however, the agreement between the pulsed and desired dynamics is very good for relatively long 
times, giving a proof of principle that indeed we may manipulate the dynamics of quantum many-body systems with 
dynamical decoupling methods [57| . 

The agreement between Uw{t) and U{t) is best at the completion of cycles, as shown in the left panel and also as 
verified, but not shown, for the system in the right panel. Interestingly, however, for the system in the middle panel 
this does not always hold. In the middle panel we compare the gapless system Hnin with A = 1/2 and the gapped 
system Hyj = Hnn /A with A = 2. The applied pulses involve a single direction z and the cycle is completed after 
only two pulses [cf. Eq. In the figure, we show with circles the values of after every pulse, so each value of the 
pair Ti_2 has associated with it two curves. For short times, the performance of the sequence is best at the completion 
of the cycles, as expected, but at longer times, the scenario is reversed and the best performance appears in between 
cycles. The origin of this shift of roles, which occurs faster for longer cycle times, deserves further investigation. 
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